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Abstract 



Recently O'Connell introduced an interacting diffusive particle system in order to 
study a directed polymer model in 1+1 dimensions. The infinitesimal generator of 
the process is a harmonic transform of the quantum Toda-lattice Hamiltonian by the 
Whittaker function. As a physical interpretation of this construction, we show that the 
O'Connell process without drift is realized as a system of mutually killing Brownian 
motions conditioned that all particles survive forever. When the characteristic length of 
interaction killing other particles goes to zero, the process is reduced to the noncolliding 
Brownian motion (the Dyson model). 

Keywords Mutually killing Brownian motions • Survival probability • Quantum Toda 
lattice • Whittaker functions • The Dyson model 

1 Introduction 

In this paper we introduce a system of finite number of one- dimensional Brownian motions 
which kill each other, and evaluate long-term asymptotics of the probability that all particles 
survive. Then we define a process of mutually killing Brownian motions conditioned that all 
particles survive forever. We show that this conditional process is equivalent to a special 
case of the process recently introduced by O'Connell in order to analyze a directed polymer 
model in 1+1 dimensions [28]. 

As an introduction of our study of many-particle systems, here we consider a family of 
one-particle systems with a parameter ,^ > 0. Let B{t),t > be such a one-dimensional 
Brownian motion that its survival probability P{t) decays following the equation 
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dP{t) 



V{B{t))Pit), t>0, 




dt 
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with a decay rate function 

^eR. (1.2) 

The function (1.2) imphes that, if the Brownian particle moves in the positive region far 
from the origin a; ^ decay of survival probability is negligible, while as it approaches 
the origin the decay becomes large. Note that the Brownian particle is able to penetrate a 
negative region x < 0, but there the decay rate of survival probability grows exponentially 
as a function of \x\. The parameter ^ is the characteristic length of the interaction to kill 
a particle acting from the origin and it is also the penetration length of a particle into the 
negative region. See Fig. 1(a). For < t < oo, if a path of the Brownian motion up to time 
t is given as {B{s) ■ < s <t}, the survival probability of the particle at the time t is given 
by an integration of (1.1), 

Pit\{B{s) : < s < 0) = exp | - ^ V{B{s))cls^ 

- exp{-^£e-^(^)/^d.}. (1.3) 

Provided that B(0) > 0, in the limit ^ — ?> 0, the process becomes an absorbing Brownian 
motion in the positive region with an absorbing wall at the origin, in which the survival 
probability (1.3) is zero if the particle hits the origin at any time s < t, and it is one if it 
stays in the positive region in the time period [0,t]. See Fig. 1(b). 

Let Q{t, y\x) be the transition probability density of this killing Brownian motion from a 
position a; e R to a position y G M during time t > 0. It will be obtained by taking an average 
of (1.3) over all realizations of Brownian path under the conditions -B(O) = y,B{t) = x. 
(Note that the time ordering is reversed, since the equation of Q given by (1.5) below is the 
backward Kolmogorov equation.) That is. 



Q{t,y\x)^E 



exp I - 2* V{B{s))ds^ 1(5(0) = y, B{t) = x) 



(1.4) 



where E[ • ] denotes an expectation with respect to Brownian motion, and 1(0;) is the 
indicator function of a condition w; l(w) = 1 if w is satisfied, and l(ci;) = otherwise. It is 
the Feynman-Kac formula (see, e.g., [15]) for the unique solution of the diffusion equation 

|g(t, y\x)={^^- V{x)^ Q{t, y\x) (1.5) 

under the initial condition Q{0, y\x) = 6{x — y). 

Let Ii,{z), e C be the modified Bessel function of the first kind 



k=Q 



^-(^)-E r(fe + i)r(fe + . + i) ' N<°°' l-^g^K-' (1-6) 
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(a) 




(b) 



Figure 1: (a) An illustration of a path of surviving Brownian particle in one- dimension with 
a killing term —V{x) given by (1.2). The parameter ^ is the characteristic length of the 
interaction to kill a particle acting from the origin. It also indicates the penetration length 
of a particle into the negative region, (b) An illustration of a path of surviving Brownian 
particle in the absorbing Brownian motion with an absorbing wall at the origin. It is regarded 
as the limit ^ — )■ of the system shown in (a). 



and Ki,{z),z G C be Macdonald's function [37, 25] defined by 

Til-Jz) - iJz) 



Kjz) 



and for integers u = n, 



sm z/vr 



|argz| < TT, for z/ 7^ 0, ±1, ±2, . . . 



KJz) 



limK^iz), n = 0,±l,±2, 



[1.7] 



;i.8) 



which are both analytic functions of z for all z in the complex plane C cut along the negative 
real axis, and entire functions of u. We see that Ii,{z) and K^{z) are linearly independent 
solutions of the differential equation 



d'^w 1 dw 



dz'^ z dz 



1 + 



w = 0. 



For X > and z/ > 0, Iy{x) is a positive function which increases monotonically as x — )■ 00, 
while K^{x) is a positive function which decreases monotonically as x 00. Then an 
integral representation of Q{t^y\x) is given by [26, 27, 28, 17] 



TT 



:i.9) 



where i = 1/— 1 and T{z) is the gamma function. 
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For each initial position a; e R, the survival probability of this killing Brownian motion 
at time T > is obtained by integrating Q{T, y\x) over y G R, 

/oo 
Q{T,y\x)dy, x eW, T > 0. (1.10) 
-oo 

If we condition the process to survive up to time T > 0, the transition probability density 
from the state (s, x) to {t, y) of the killing Brownian motion is given by 

P^{s,x;t,y) = ^^,^^Q{t-s,y\x), x,y eR, 0<s<t<T. (1.11) 

If we use the integral representations of Ki,{z) [37, 25], 

, , 2'^r(z/ + l/2) cos(zu) , 

and 

K.{z) = ]^£ s^-' exp + 1^ I ds, (1.13) 

we can evaluate the long-term asymptotics of the survival probability for \x\ < oo, 

A^(T,x)=Cr-V2Xo(e-^/^)x|l + o(^-^)} in (1.14) 

with C = 3^-^27 TT as shown in Appendix A. Then we can take the limit T — > oo of (1.11) as 

P(t,y\x) = lim P^{0,x;t,y) 

= -jP—rj^MtM^), x,yeR, 0<i<oo. (1.15) 

This is the transition probability density of the killing Brownian motion conditioned to 
survive forever. We should note that this conditional Brownian motion is the diffusion 
process with infinitesimal generator 

Matsumoto and Yor [26, 27] have studied the exponential functionals of Brownian motion 

ZuY{t) ^ ilog^^^ j\'^^'^'^ds^-B{t), t>0, (1.17) 

which is the diffusion process whose infinitesimal generator is (1.16). We can say that the 
Matsumoto- Yor process (1.17) is realized as the present killing Brownian motion conditioned 
to survive forever [17]. 
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Now we consider the limit ^ ^ of the formulas obtained above. Since Iiy{z) is defined 
as the series expansion (1.6), we can see /j^fc(e~^/^) ~ {e~^/^ /2y^^ /T{i^k + 1) ~ g-ifex 
{ ^ for a; > 0. Then (1.7) gives Ki^k{e~'''^) ^ sm{kx)/ {^k)l{x > 0) in ^ ^ 0. Therefore 
the integral of (1.9) can be performed in the limit and we have 

q{t,y\x) = \imQ{t,y\x) 

= -^|e-(--^)V2* _ e-(-+3')V2A > 0,y > 0). (1.18) 

It is the transition probability density of the absorbing Brownian motion with an absorbing 
wall at the origin, which is easily obtained by applying the reflection principle of Brownian 
motion [15]. Applying the result (A. 2) shown in Appendix A, we see that 

lim^i^o(e"'^/^) = X l(x > 0), (1.19) 
and then (1.15) gives in this limit 



pit,y\x) = lim P{t,y\x) 

= -q{t,y\x), x>0,y>0, t>0, (1.20) 

X 

which is a harmonic transform (/i-transform) of (1-18) and is identified with the transition 
probability density of the three-dimensional Bessel process, BES(3) (see, for instance, [20]). 
As a matter of fact [26, 27], (1.17) gives lim^^o ZuYit) = 2 supQ<g<j B{s) —B{t),t > 0, which 
is indeed distributed as BES(3) (Pitman's 2M — X theorem [30]). 

In the present paper, we consider an iV-particle system of one-dimensional Brownian 
motions with N > 2, B{t) — {Bi{t), . . . , BN{t)),t > 0, with a positive parameter ^ > 0, 
such that the probability that all N particles survive up to time t, i-*iv(^), decays following 
the equation 

dPNit) 
dt 

with a decay rate function 



^-V:,{B{t))PN{t), t>0 (1.21) 



7V-1 

Vn{x) = -t^Y1 e-(^^+^-^^)/^ X e R^. (1.22) 

We study an integral representation of the transition probability density QN{t, ul^c), x,y E 
]R^,t > 0, which is a multivariate extension of (1.9). The extensions of the formulas (1.14) 
and (1.15) are shown. We prove that the A^- particle system of the mutually killing Brow- 
nian motions following (1.21) and (1.22) conditioned that all N particles survive forever is 
equivalent to a special case (without drift, i/ = 0) of the O'Connell process [28, 17]. 

In the context of quantum mechanics, the decay rate functions of survival probability 
(1.2) and (1-22) in the killing Brownian motions are considered to giving potential energy 
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of the systems. They are called the Yukawa potential and the Toda-lattice potential [36], 
respectively. The multivariate extensions of Macdonald's function (1.7) with (1.8) are the 
Whittakcr functions (see [1, 28, 4, 2, 29] and references therein), which have been extensively 
studied in mathematical physics as eigenfunctions of the quantum Toda lattice [24, 32, 31, 
9, 21, 22, 23, 14, 7, 8]. See also [10]. 

In Sect. 2, as preliminaries, useful integral representations of the Whittaker functions are 
given for the GL{N, R)-quantum Toda lattice. The transition probability density of the 
system of mutually killing Brownian motions is given, in the case that killing of particles 
does not occur during an observing time period, as an integral of a product of the Whittaker 
functions over the Sklyanin measure. Then the main results are given. We also demonstrate 
that when the characteristic length ^ of the interaction kilhng other particles goes to zero, 
the O'Connell process is reduced to the noncoUiding Brownian motion, which is known to be 
equivalent to the Dyson model [11, 20] for the eigenvalue process of Hermitian matrix- valued 
diffusion process {i.e. Dyson's Brownian motion model with the parameter /3 = 2 [5]). The 
proofs of Lemmas and Proposition are given in Sect. 3. Appendix is given for proving (1.14) 
and (1.19) used above. 



2 Preliminaries and Main Results 

2.1 Quantum Toda lattice and Whittaker functions 

In this subsection we set ^ = 1 in Eq.(1.22) and consider Vn{x) as a potential energy of 
a quantum A^-particle system in one dimension. Then we have the Hamiltonian of the 
GL{N, R)-quantum Toda lattice 

1 ^ pp. 

^^ = -2Ea^ + E^"^^^"^~^^^ (2-1) 

j=l j=l 

for e {2, 3, . . . }. (In this paper, we set the mass of particle m — 1 and H—IIot: simplicity 
of notation.) The generalized eigenvalue problem 

HN'^xix) = \^x{x), X e M^, (2.2) 

is solved by setting 

1 ^ 1 

A = --$]^| = -2M^ u^{ly^,...,u^)eC^ (2.3) 



with the eigenfunction "ifxix) — ipi, (x), where ipl, '{x) is the GL{N, R)-Whittaker function 
[24, 32, 31, 21] (the class-one gt^- Whittaker function [8, 1, 28, 4, 2]). When = 2, it is 
expressed by using Macdonald's function (1.7) as 

^S,,,)(xi,a;2) = 2e(^i+'^^)(^^+^^)/^X,,_,,(2e-(^-^i)/^). 
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In this sense, the Whittaker functions {'ipiy \x)} jsf>2 are regarded as multivariate extensions 
of Macdonald's function. 

Corresponding to the two kinds of integral representations (1.12) and (1.13) for K^{z), 
two integral representations arc known for the Whittaker function 'iIj^j^\x),u G C^. The 
integral expression corresponding to (1.12) is the classical one originally given by Jacquct [13] 
and was rewritten as the following form by Stade [33]. (Here we use the notation given by 
[21].) Let Z denote an upper triangular N xN matrix with unit diagonal; Z = (^j,fc)i<j,fc<Af, 



1, 

0, 



l<j<k<N, 
l<j = k<N, 
l<k<j<N. 



(2.4) 



We write the integral of a function f of Z over all real Z as 

N k-l 



r ^ ^ poo 

/ /(z)dz^nn/ dz,,kf{z). 



k=l j = l 



The transpose of Z is denoted by *Z' and the principal minor of size j of matrix Z^Z is 
written as Aj(Z*Z), 1 < j < iV. For a; = (xi, . . . , xjv) G M^, i/ = (z^i, . . . , i/jv) G C^, 

l<j<k<N ^ 

J]RiV(iV-l)/2 -^J^ I ' J 



dZ. (2.5) 



Another integral representation, which corresponds to (1.13), was given by Givental [9]. See 
also [23, 7]. Let T denote a lower triangular array with size N, T = {Tj^k, ^ k < j < N). 
For a given x e M^, let Tn{x) be the space of all real triangular arrays T with size N 
conditioned 

TN,k ^Xk, l<k<N. (2.6) 
We write the integral of a function / of T over T^ix) as 

/ f{T)dT^llll dT,,,fiT)ll6iT^,e-xe). 



Then 



exp 



TNiX) 

N-1 j 



N 



i-1 



j=l \k=l k=l 

— ^^-('^j,k-Tj+i,k) _|_ g-(rj+i_fc+i-rj_fc) I 



j=i k=i 



dT. 



(2.7) 
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Remcirk 1. There is the 'third' kind of integral representation of the Whittaker function, 



N-1 

n 



ui=iUSniw-i,k) 



N 



X 



^ N-1 j 

exp <j i ^ J2i'yj,k - lj-i,k)xj \ n n 

j=l k=l ) j=l k=l 



(2.8) 



where (7Ar,i, . . . , 7Af,Ar) = (z^i, . . . , i^n), lj,k = for A; > j, and the domain of integration S is 
defined by the conditions minfc{Im7j^fc} > maxfc{Im7j+i^/;} for aU j = 1, 2, . . . , — 1. This is 
called the Mellin-Barnes integral representation, since it can be regarded as the multivariate 
extension of the MeUin-Barnes representation of Macdonald's function [37] 



where the path of integration C being a vertical line to the right of any poles of the integrand. 
The derivation of the MeUin-Barnes integral representation (2.8) from the classical one (2.5) 
is shown in [34, 12] (see also [22, 23]). The equivalence between the Givental integral rep- 
resentation (2.7) and the MeUin-Barnes integral representation (2.8) is fully discussed in 
[8]. 



2.2 Transition probability density of the mutually killing Brown- 
ian motions 

For e {2, 3, . . . }, we consider the A^-particle system of mutually killing Brownian motions, 
in which the probability that all particles survive decays in time following (1.21) with 
(1.22) depending on realization of paths. For x,y G M^,t > 0, the transition probability 
density from the state x to y during time interval t is then given by 



QN{t,y\x) = E 



E 



exp j- ^ VNiB{s))dsj 1(B(0) = y, B{t) = x) 

^*e-(^^+i(^)-^^(^»/«ds| 1(5(0) = y, B{t) = x) 



■ (2.9) 



It is the Feynman-Kac formula for the solution of the diffusion equation (the backward 
Kolmogorov equation) 

= ^Au{t, x) - V^{x)u{t, x) (2.10) 



with A = d'^/dx^ under the initial condition 



N 



U 



(0, x) ^5{x-y) = \[ 5{xj - Vj). 



(2.11) 
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X 



X 



A 



(b) 



Figure 2: (1) An illustration of three paths of surviving Brownian particles in the mutually 
killing Brownian motions. The effect of killing is not negligible when the distance between the 
nearest neighbor pair of particles is less than ^; Bj^iit) — Bj(t) < ^. The ordering of particle 
positions can be changed, but the risk of pair annihilation becomes large exponentially, when 
Bj(t) — Bj^iit) > ^. (b) In the hmit ^ — > 0, the present system converges to the vicious 
Brownian motion, in which the process survives if and only if paths are not intersecting. 

Remark 2. If we regard (2.10) as a diffusion equation describing an A^- dimensional Brow- 
nian motion in with a killing term —Vn{x), (2.9) gives a transition probability density 
for a particle assumed to be at a position x G such that it survives during time t and 
it arrives at a position y G after the time duration t (the Feynman-Kac formula, see, 
for instance, [15]). In the present paper, on the other hand, we would like to consider an 
A^-particle system of one-dimensional Brownian motions, such that the j-th and (j + l)-th 
particles will be pair annihilated with high probability if Bj^i{t) — Bj{t) < ^,1 < j < N — 1, 
and then the probability that all particles survive avoiding from any mutual killing de- 
cays in time following (1.21) with (1.22). In order to discuss processes, in which mutual 
killing of particles actually occurs and total number of particles decreases in time, we have 
to specify the stochastic rule to determine which pair of particles is annihilated, when the 
survival probability conditional on a path, exp{— VN{B{s))ds}, becomes small. Here we 
are interested in, however, the situation that mutual killing of particles does not occur at 
all and all A^ particles survive, following the notion of vicious walker models [6, 18, 3, 17]. 
The Feynman-Kac formula (2.9) gives the transition probability density between A^-particle 
configurations x,y & M^. Figure 2(a) illustrates three paths of surviving particles, in which 
change of ordering of particle positions occurs within the spatial scale ~ ^. 

By the fact that the Whittaker function iIj^^\x) solves (2.2) with (2.3), we can see that 




N 



(2.12) 
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solves the diffusion equation (2.10). 

Tfie density function of the Sklyanin measure [32] is defined by 



"^(^^ = 72^ n \mh-m-' (2.13) 

^ ' ' i<j<e<N 



for k e C^. By Euler's reflection formula r(z)r(l - z) ^ 7r/sin(7rz), if fc e , 



i<j<e<N i<j<e<N 
Since the Whittaker functions satisfy the completeness relation with respect to the Sklyanin 



measure SN{k)dk = SN{k)Y[f=idkj [31, 21, 23, 8], 



/ ^l;f{x)^'-''Uy)s^{k)dk^6{x-y), a;,yeM^ (2.14) 
the integral of the solution (2.12) multiplied by V^'^^^^^dz/O on this measure, 

Qr,{t, y\x) = e-'\f'\'/'^l^jl{x/0^^^lf.{y/OsNmdk (2.15) 

satisfies the initial condition (2.11). That is, (2.15) is an integral representation of the 
transition probabiUty density (2.9). It should be noted that the orthogonality relation [21, 8] 

/ <fc(^)^^fc'(^)^^ = ^^E^(^-^(^'))' (2.16) 

JjRAT -^K SN{k)N\ ^ 

is established, where Sn is a set of all permutation of N indices and cr(fe') = (^^(i), • • • , ^^(Af))- 
It guarantees the Chapman-Kolmogorov equation for the transition probability density 

/ QN{t2.z\ii)Q^{t^,y\x)dy^Q^{tr^t^,z\x), a;,zeR^ (2.17) 

for < ti, ^2 < oo, which should hold, since the present process is Markovian. 

Remark 3. The completeness relation (2.14) of the Whittaker functions is a multivariate 
extension of the relation for Macdonald's functions, 

^ [ Ki„{e-^)K_ik{e-y)ksmh{7ik)dk = d{x-y), x,y eR, (2.18) 

which will define the Kontorovich-Lebedev transform [25]. We can use (2.18) to confirm that 
(1.9) satisfies the initial condition Q{0,y\x) — S{x — y). 
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2.3 Asymptotics of transition probability density and survival prob- 
ability 

By using the integral representation (2.5), we can prove the following asymptotics of the 
transition probability density. 

Lemma 2.1 For x,y e , \x\ < oo,t > 0, as \x\/y/i 0, 

Af2/2 



Jr^ i<j<ifc<Ar 



\y\y2t 



\x\ 



On the other hand, by using Givental's integral representation (2.7), we have the following 
[28]. Let Wjv be the Weyl chamber of type Ajv-i, 

Wat — {x — (xi, X2,..., Xn) £ : Xi < X2 < ■ ■ ■ < Xn} 

and set 

hN{x)= Yl (Xk-Xj). (2.20) 

l<j<k<N 



Lemma 2.2 For x e 



hm /3-^(^-)/Vr(/3a.) = ^^Ux G W^). (2.21) 

n,=i j! 

The probability that all N particles survive up to time t in this system of the mutually 
killing Brownian motions is given by 

J\fN{t,x)^ I QN{t,y\x)dy (2.22) 

for each initial state x G M^. We call it simply the survival probability in this paper. The 
following long-term asymptotics of the survival probability is obtained. 

Proposition 2.3 Let 

Cn = e"^"-^^/^ ' . A, (2.23) 

7^ 11^=1 j! 

with 

A^^j dri f d/ue-d^l^+l^l^) n \iVk - V^ipik - Pij)' + iVk - V^f}]. (2.24) 

^W^ JWm l<j<k<N ^ 

Then for xeR^, 

lim i^(^-^)/Wiv(t, x) = Cn4''\^/0- (2-25) 

t— >-oo 
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2.4 O'Connell process as a system of mutually killing Brownian 
motions conditioned that all particles survive forever 

Let < T < oo and we consider the A^-particle system of the mutually killing Brownian 
motions conditioned that all particles survive up to time T. For < s <t <T, the transition 
probability density from the state {s,x) to {t,y),x,y e of this conditional process is 
given by 

Pl{s, X- 1, y) = ^I^ll ~ ^' ^\ QN{t - s, y\x), 0<s<t<T. (2.26) 

JMNyl — S, X) 

The fact that (2.26) depends not only t — s but also T — t and T — t implies that this 
conditional process is inhomogeneous in time. 

Proposition 2.3 enables us to take the hmit T — > oo of (2.26). As a result we obtain the 
temporally homogeneous process with the transition probability density 

PN{t,y\x) = lim Pj^{0,x;t,y) 

T^oo 

= ^^^Dr^QN{t,y\x), x,yeR\t>0. (2.27) 

It is easy to confirm the following [17]. 

Proposition 2.4 The transition probability density (2.27) of the N -particle system, of the 
mutually killing Brownian motions conditioned that all particles survive forever solves the 
following diffusion equation, 

^u(t, x) = ^Aii(i, x) + V log Vr)(a;/0 • Vu{t, x) 

= iE^''('.-)+ E '-/^' ^,.(,.). .eR~.*>0. (2.28) 

j=i i i<j<N ^ ^ 



under the initial condition u{0, x) — S{x — y),y & 



Recently O'Connell [28] introduced a diffusion process in with the infinitesimal gen- 
erator of the process 



- -(#)-^(^^ + IH>r 



= ^A + VlogV'if^(a;)-V (2.29) 
with drift u e M^. Then we can state the following. 

Theorem 2.5 The O'Connell process without drift, u — 0, is realized as the system of 
mutually killing Brownian motions conditioned that all particles survive forever. 
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In an earlier paper [17], we considered the system of mutually killing Brownian motions with 
drifts = (i/i, . . . , un) e conditioned that all particles survive forever and took the hmit 
1/ ^ to realize the O'Conncll process without drift i/ = 0. In the present paper we have 
shown that an introduction of drift u is not necessary to deriving (2.27) and concluding 
Theorem 2.5. 

2.5 The ^ ^ limit 

It is known that the functions {m'^^^x, k)}f^^^i^ called the fundamental Whittakcr functions 
are defined so that the (class-one) Whittaker functions discussed so far are expressed by the 
following alternating sum of them [21, 1], 

l<j<£<N ^ * aeSN 

where (j{k) — (/co-(i), ■ ■ ■ , ka{N)) for each permutation a E Sn- The functions {m^^\x, k)}j^^^N 
are normalized as 

limm^^\x/^,m = e'^-^ for x G Wat- (2.31) 

Note that (2.30) is a multivariate extension of the definition (1.7) of Ki,(z), which is given 
by an alternating sum of and I_^{z). Therefore, when x,y & Wjv, we can see that 

(2.15) with (2.13) gives 

qN(t,y\x) = \imQN{t,y\x) 

1 ^ . . /27r^ 



(2vr)^ , 



Y: sgn(p) e-l^-(2/)lV^* 



det 

l<j,fc<JV 



(2.32) 



This is the Karlin-McGrcgor determinant giving a total mass of nonintersecting paths of 
Brownian motions from x G Wjv to y E during time t [16, 20]. As a matter of fact, 
when x,y G Wat, the limit ^ — >■ of (2.9) should give the transition probability density 
of the A'"-dimensional absorbing Brownian motion in the Weyl chamber W^f (the vicious 
Brownian motion [17]). See Fig.2(b). Lemma 2.2 implies that, if a;, y G Wjv, 

PN{t,y\x) = limPN{t,y\x) 

= ^^^qN{t,y\x), t>0, (2.33) 
hN{x) 
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which is an /i-transform of qn by the Vandermonde determinant (2.20). It shows that 
the hmit ^ — > of the O'Connell process is identified with the Dyson model (Dyson's 
Brownian motion model with (3 = 2 [5]), which is the A'^-particle system of Brownian motions 
conditioned never to collide with each other (the noncoUiding Brownian motion) [18, 19, 20]. 

Remark 4. The survival probability up to time T < cxd of the one-dimensional absorbing 
Brownian motion with an absorbing wall at the origin is given by 

/oo 
q{T,y\x)dy, (2.34) 
-oo 

where the transition probability density q is given by (1.18). The long-term behavior is 
readily obtained as 

AA°(r, x) ~ c^xT-'f', x>0, T ^ oo (2.35) 

with c° = •\/2/7r, = 1/2. Therefore, the exponent is the same with that found in the 
asymptotics (1-14) of the survival probability of the killing Brownian motion discussed in 
Sect.l; 

^^{T, x) ~ 3^^Ko{e-''/^)T-"^, x eR, T ^ oo. (2.36) 
If we take the limit ^ ^ 0, (1-19) gives 

limA/'(r, x) ~ cxT-"^, x>0, T^oo (2.37) 

with c = 3-\/2/7r = 3c°. For N > 2, Proposition 2.3 imphes 

Xn{T, x) ~ Cjv4^\a;/0T-'^^, xeR^, T ^ oo (2.38) 
with Cat given by (2.23) and with the survival probability exponent 

0iv = \n{N - 1) (2.39) 

for the present A'^-particle system of mutually killing Brownian motion. This exponent is 
the same as that for the vicious Brownian motion (the absorbing Brownian motion in W^^) 
[6, 18]; 

AA°(r,a;) = / qj,{T,y\x)dy 

~ c%hNix)T-'f"' , xeWN, T^oo (2.40) 

with 



„o _ 2 ^/^nf=ir(j72 ) 
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Lemma 2.2 gives then 



limTVivfT, x) ~ CNhN(x)T-'f"' , x e Wjv, T ^ oo, (2.42) 

where 

23Af(Af-l)/4 

with (2.24). The above results (P and cat 7^ c^, N >2 are consequences of the facts that 
for the transition probabihty densities Q{t^y\x) and QN{t,y\x),N > 2, the long-term limit 
T — )■ cxD and the limit ^ ^ are noncommutable. 

3 Proofs of Lemmas and Proposition 
3.1 Proof of Lemma 2.1 

By the integral representation (2.5), we put i'^.^^{x/^) = exp(ifc ■ x)'i//^^{x / ^) . Then by 

changing the integral variables as iJ,j = ^JTj2kj — i{xj — yj)/\/2i, I < j < N, (2.15) with 
(2.13) is written as 

-2 

l<j<k<N 



^ i^e-^^^'^iaUi^/Oi^%t^)iy/0 n |r(i(«.(/..) - «,(/.,))) 

where 



and Q!(/^) = (ai(^i), . . . , Q;Ar(/iAr)). For \x\ < 00, in the limit \x\/\^ — )■ 0, aji^J^j) — 
iVVt{H - iyj/V2i} ^ 0, 1 < J < AT, and thus ^£}^)(x/C) ^ ipi^\x/0 = ^^^^a^/O- 
Since 2;r(2;) — >■ 1 as 2; — > 0, 

-2 2 f 1 1 

r(i(Q;fc(///c) - oijinj))) ~ J _ ^^.^2 + — (y^ - y^.)2 1 

in this limit. Then (2.19) is obtained. | 
3.2 Proof of Lemma 2.2 

The following argument is found in Section 6 of O'Connell [28]. For N > 2, the Gelfand- 
Tsetlin pattern {complete branching) is defined as [35] 

GTn = |(a;^a;^...,a;^) e Wi X W2 X ••• X W;v : 

4"' < 4 < 4X1 l<k<j<N-l], (3.1) 
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where Wj = {x = {xi, . . . ,Xj) eW : xi < ■ ■ ■ < xj}, I < j < N. By the Givental integral 
representation (2.7), 



v^r(/3a.) = / exp 



N-l j 



{Tj,k-Tj+i,k) _|_ g-(^3+i,fc+i-?), 



^iV(iV-l)/2 

'r;v(a;) 



N-l j 



|g-/3(T"'j,fe-T"j+i,fe) _|_ g-/3(r'j+i^fe+i-T'j^i.) I 



j = l fc=l 



dT', 



where we set T ^ ^T' ^ {l3T'j,k ■l<k<j<N). Let t^' = (T'^- 1, . . . , T'jj), l<j<N.As 
/3 ^ oo, if 

e GT^ and t^ = xeWN, (3.2) 

then the integrand converges to 1, and otherwise it becomes zero. Then the integral is the 
volume of GT^v under the condition — x & Wn (the GT-polytope), which is given by 
hN{x)/ Y[f=i j^-- Then the proof is completed. | 

3.3 Proof of Proposition 2.3 

By definition (2.22), Lemma 2.1 gives 



J\fN{t,x) = C 



N{N-1)_ 



1 



l<j<k<N 

in \x\/\/i — )■ 0. If we put yj/\/2i = r]j, 1 < j < N, it is written as 



\x\ 



MNit.x) = C 



N{N-1)_ 



1 



m V t 



N{N~l)/2 



X 



X 



l<j<k<N 

in \x\j\/t 0. By Lemma 2.2, 



n [(/^.-/^.)' + (%-r^,)'] x|l + o(^M^|, 



( 



e 7 V ^ 

in f — > oo, and then the proposition is proved. 



N{N-l)/2 

priv 1 / (^ e Wat), 

lii=l 
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Appendix 



A One-dimensional killing Brownian motion 

By (1.12), we can put ^^^^(e-^/^) = e^'=^X,^fe(e-^/^) with K,{z) = {2''r{iy + 1/2)/^) 
du cos{zu) / {1 + 14^)'^+^/^. Changing the integral variable in (1.9) by — i{x — 

y) I — fj,, the transition probability density can be written as 

71 \ t 

e-^'i?ia(;.) (e-^/«)i?_,«(^) (e-^/«) |r(za(/.)) 



with Q;(yu) = ^y^2/t{fj, + i{x — y)l\f2±\. For < oo, as x/^/t ^ 0, q;(/x) — ^-^2/t(/i — 
iyl\f2t) 0, and r(m(/x)) ~ {ii^l2ft{ix - iyl\f2t)Y^. Then 



g(t, yk) = ^ ( - J e-^'/2*/ro(e-^/«)i^o(e-^/«) 
= e^t-'/'e-y'/'^ (^1 + ^) Ko(e"-/«)i^o(e 



in a;/\/t ^ 0, since /To = Kq by definition. It gives through (1.10) 



(A.l) 



with 



= J\-y'/''(^l + y-PjKo{e-y/^)dy 

/oo 
e-"'(l + 2M2)i^o(e-^"/«)dM. 
-oo 

By (1.13), for /3 > 0, /s:o(e-'^'=) is written as 
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where we have set s — e^'". In the hmit ^ — >■ oo, e + e ^(^+^) — 0, ii v < x and v > —x, 

and e"^^^"""^ + e"^^^'^'"^ — oo, otherwise. Therefore 

hm /3-^Ko{e-^'') = - / dvl{x >0)^xl{x>0). (A.2) 

Then i^o(e"^"/^) V2tM/Cl(M > 0) in t ^ oo, and 

It = -rl e-''\l + 2u^)udux{l + o{l)} = — x{l + o{l)} in t ^ oo. 
Put this result into (A.l), then we obtain the asymptotics (1.14). 

Acknowledgements The present author would like to thank T. Sasamoto and T. Ima- 
mura for useful discussion on the present work. A part of the present work was done during 
the participation of the present author in Ecole de Physique des Houches on " Vicious Walkers 
and Random Matrices" (May 16-27, 2011). The author thanks G. Schehr, C. Donati-Martin, 
and S. Peche for well-organization of the school. This work is supported in part by the 
Grant-in-Aid for Scientific Research (G) (No. 21540397) of Japan Society for the Promotion 
of Science. 

References 

[1] Baudoin, P., O'Gonnell, N.: Exponential functionals of Brownian motion and class-one 
Whittaker functions. Ann. Inst. H. Poincare, B 47, 1096-1120 (2011) 

[2] Borodin, A., Gorwin, I.: Macdonald processes. arXiv:math. PR/1111. 4408 

[3] Gardy, J., Katori, M.: Families of vicious walkers. J. Phys. A 36, 609-629 (2003) 

[4] Gorwin, L, O'Gonnell, N., Seppalainen, T., Zygouras, N.: Tropical combinatorics and 
Whittaker functions. arXiv:math. PR/1110.3489 

[5] Dyson, F. J.: A Brownian-motion model for the eigenvalues of a random matrix. J. 
Math. Phys. 3, 1191-1198 (1962) 

[6] Fisher, M. E.: Walks, walls, wetting, and melting. J. Stat. Phys. 34, 667-729 (1984) 

[7] Gerasimov, A., Kharchev, S., Lebedev, D., Oblezin, S.: On a Gauss-Givental represen- 
tations of quantum Toda chain wave equation. Int. Math. Res. Not. 1-23 (2006) 

[8] Gerasimov, A., Lebedev, D., Oblezin, S.: Baxter operator and Archimedean Hecke 
algebra. Gommun. Math. Phys. 284, 867-896 (2008) 

[9] Givental, A.: Stationary phase integrals, quantum Toda lattices, flag manifolds and 
the mirror conjecture. In: Topics in Singular Theory, AMS Trans. Ser. 2, vol. 180, 
pp.103-115, AMS, Rhode Island (1997) 



18 



[10] Gorsky, A., Nechaev, S., Santachiara, R., Schehr, G.: Random ballistic growth and 
diffusion in symmetric spaces. arXiv:math-ph/1110.3524 

[11] Grabiner, D. J.: Brownian motion in a Weyl chamber, non-colliding particles, and 
random matrices, Ann. Inst. Henri Poincare, Probab. Stat. 35, 177-204 (1999) 

[12] Ishii, T., Stade, E.: New formulas for Whittaker functions on GL(n, M). J. Func. Anal. 
244, 289-314 (2007) 

[13] Jacquet, H.: Fonctions de Whittaker associees aux groupes de Chevalley. Bull. Soc. 
Math. France 95, 243-309 (1967) 

[14] Joe, D., Kim, B.: Equivariant mirrors and the Virasoro conjecture for flag manifolds. 
Int. Math. Res. Notices, no.l5, 859-882 (2003) 

[15] Karatzas, I., Shreve, S. E.: Brownian Motion and Stochastic Calculus. 2nd edn. 
Springer, (1991) 

[16] Karlin, S., McGregor, J.: Coincidence probabilities. Pacific J. Math. 9, 1141-1164 (1959) 

[17] Katori, M.: O'Connell's process as a vicious Brownian motion, to appear in Phys. Rev. 
E. arXiv:math-ph/1110.1845 

[18] Katori, M., Tancmura, H.: Scaling limit of vicious walks and two-matrix model. Phys. 
Rev. E 66, 011105 (2002) 

[19] Katori, M., Tanemura, H.: Non-equilibrium dynamics of Dyson's model with an infinite 
number of particles. Commun. Math. Phys. 293, 469-497 (2010) 

[20] Katori, M., Tanemura, H.: NoncoUiding processes, matrix-valued processes and deter- 
minantal processes. Sugaku Expositions 24, 263-289 (2011); arXiv:math. PR/1005.0533 

[21] Kharchev, S., Lebedev, D.: Integral representation for the eigenfunctions of a quantum 
periodic Toda chain. Lett. Math. Phys. 50, 55-77 (1999) 

[22] Kharchev, S., Lebedev, D.: Eigenfunctions of GL(iV, M) Toda chain: the Mellin-Barnes 
representation. JETP Lett 71, 235-238 (2000) 

[23] Kharchev, S., Lebedev, D.: Integral representations for the eigenfunctions of quantum 
open and periodic Toda chains from the QISM formalism. J. Phys. A: Math. Gen. 34, 
2247-2258 (2001) 

[24] Kostant, B.: Quantisation and representation theory. In: Representation Theory of Lie 
Groups, Proc. SRC/LMS Research Symposium, Oxford 1977, LMS Lecture Notes 34, 
pp. 287-316, Cambridge University Press, Cambridge (1977) 

[25] Lebedev, N. N.: Special Functions and Their Applications. Prentice-Hall, Inc. (1965) 



19 



[26] Matsumoto, H., Yor, M.: An analogue of Pitman's 2M — X theorem for exponential 
Wiener functional, Part I: A time-inversion approach. Nagoya Math. J. 159, 125-166 
(2000) 

[27] Matsumoto, H., Yor, M.: Exponential functionals of Brownian motion. I: Probability 
laws at fixed time. Probab. Surveys 2, 312-347 (2005) 

[28] O'Connell, N.: Directed polymers and the quantum Toda lattice, to appear in Ann. 
Probab. arXivimath. PR/0910.0069 

[29] O'Connell, N.: Whittaker functions and related stochastic processes. 
arXiv:math. PR/1201.4849 

[30] Pitman, J. W.: Onc-dimensional Brownian motion and the three-dimensional Bessel 
process. Adv. Appl. Prob. 7, 511-526 (1975) 

[31] Semenov-Tian-Shansky, M.: Quantisation of open Toda lattices. In: Dynamical Systems 
Vll : Intcgrable Systems, Nonholonomic Dynamical Systems. Edited by V. 1. Arnol'd 
and S. P. Novikov. Encyclopedia of Mathematical Sciences, vol.16. Springer, Berlin 
(1994). 

[32] Sklyanin, E. K.: The quantum Toda chain. In: Non-linear Equations in Classical and 
Quantum Field Theory, Lect. Notes in Physics, 226, pp. 195-233, Springer, Berlin 

(1985) 

[33] Stadc, E.: On explicit integral formulas for GL(n, M)- Whittaker functions. Duke Math. 
J. 60, 313-362 (1990) 

[34] Stadc, E.: McUin transforms of GL{n,R) Whittaker functions. Amer. J. Math. 123, 
121-161 (2001) 

[35] Stanley, R. P.: Enumerative Combinatorics, vol.2, Cambridge University Press, Cam- 
bridge (1999) 

[36] Toda, M.: Theory of Nonlinear Lattices. 2nd edn. Springer, Berhn (1989) 

[37] Watson, G. N.: A Treatise on the Theory of Bessel Functions. 2nd edn. Cambridge 
University Press, Cambridge (1944) 



20 



